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ORIGINAL ARTICLE 

PKPD Modeling of VEGF, sVEGFR-2, sVEGFR-3, and sKIT 
as Predictors of Tumor Dynamics and Overall Survival 
Following Sunitinib Treatment in GIST 

EK Hansson 1 , MA Amantea 2 , P Westwood 1 , PA Milligan 2 , BE Houk 2 , J French 23 , MO Karlsson 1 and LE Friberg 1 

The predictive value of longitudinal biomarker data (vascular endothelial growth factor (VEGF), soluble VEGF receptor 
(sVEGFR)-2, sVEGFR-3, and soluble stem cell factor receptor (sKIT)) for tumor response and survival was assessed based 
on data from 303 patients with imatinib-resistant gastrointestinal stromal tumors (GIST) receiving sunitinib and/or placebo 
treatment. The longitudinal tumor size data were well characterized by a tumor growth inhibition model, which included, as 
significant descriptors of tumor size change, the model-predicted relative changes from baseline over time for sKIT (most 
significant) and sVEGFR-3, in addition to sunitinib exposure. Survival time was best described by a parametric time-to-event 
model with baseline tumor size and relative change in sVEGFR-3 over time as predictive factors. Based on the proposed 
modeling framework to link longitudinal biomarker data with overall survival using pharmacokinetic-pharmacodynamic models, 
sVEGFR-3 demonstrated the greatest predictive potential for overall survival following sunitinib treatment in GIST. 
CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e84; doi:1 0. 1 038/psp.201 3.61 ; published online 20 November 2013 



Sunitinib malate (SUTENT, Pfizer, NY) is an oral multitargeted 
tyrosine kinase inhibitor with antitumor and antiangiogenic 
properties. Sunitinib acts by inhibition of vascular endothelial 
growth factor receptors (VEGFR-1, VEGFR-2, and VEGFR- 
3), platelet-derived growth factor receptors (PDGFR-a and 
PDGFR-J3), stem cell factor receptor (KIT), Fms-like tyrosine 
kinase-3 receptor (FLT3), colony stimulating factor receptor 
type 1 (CSF-1 R), and the glial cell line derived neurotrophic 
factor receptor (RET). 1-4 Sunitinib is currently approved 
in e.g., the United States and Europe for the treatment of 
advanced renal cell carcinoma, imatinib-resistant gastroin- 
testinal stromal tumors (GIST), and pancreatic neuroendo- 
crine tumors. 

With the introduction of new approaches for treating can- 
cer, i.e., targeted therapies, the need for new approaches to 
identify an effective dose and to assess treatment response 
has evolved. Using the traditional concept of maximum-toler- 
ated dose is not appropriate for these agents because clini- 
cally effective dose could be achieved in advance of manifest 
toxicity. 5 Furthermore, the use of the standard Response 
Evaluation Criteria in Solid Tumors 6 and tumor shrinkage as 
criteria for treatment response may be problematic. Antian- 
giogenic drugs are typically cytostatic in their mechanism of 
action and consequently tumor reduction may be a less valu- 
able index of efficacy. 5 

Identification of (causal path) biomarkers could enable 
dose optimization and monitoring of response based on 
the changes in levels. The identified biomarker relationships 
could also be used to increase the understanding of the 
mechanism of action, demonstrate proof of concept in early 
phases of drug development, and enable individualization of 
ongoing treatment. 7 



Mechanism-based biomarkers known to change with 
sunitinib treatment are the vascular endothelial growth fac- 
tor (VEGF) and soluble fragments, generated by proteolytic 
cleavage, of the KIT (sKIT) and VEGF receptors (sVEGFR-2 
and sVEGFR-3). 8-11 Although tumor size at the start of treat- 
ment and changes in tumor size from baseline at weeks 7 or 
8 have previously been proposed as predictors of survival in 
various tumor types and treatments, 12-16 there is the potential 
that these tumor measures could serve as mechanism-inde- 
pendent biomarkers in sunitinib-treated patients. However, 
as stated previously, sunitinib has a cytostatic mechanism of 
action, and potentially a lower effect on tumor size could be 
expected. 17 

To identify clinically relevant outcome predictors and deter- 
mine their optimal time for measurement, an understanding of 
the underlying exposure-biomarker-effect relationship over 
time is essential. The use of pharmacokinetic-pharmacody- 
namic (PKPD) models can enable simultaneous analysis of 
both longitudinal biomarker and survival data. The entire time 
course of individual changes can then be used to assess 
relationships among biomarkers that may be separated in 
time. 18 In this analysis, exposure-effect relationships were 
characterized using nonlinear mixed-effects PKPD models to 
evaluate VEGF, sVEGFR-2, sVEGFR-3, and sKIT as poten- 
tial predictors of tumor response and subsequent overall sur- 
vival following sunitinib treatment in GIST (Figure 1). 

RESULTS 

Data on sunitinib exposure (area under the concentration- 
time curve, AUC), biomarkers, tumor size (sum of longest 
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Figure 1 Investigated relationships for the evaluation of vascular 
endothelial growth factor (VEGF), soluble VEGF receptor 
(sVEGFR)-2, sVEGFR-3, and soluble stem cell factor receptor 
(sKIT) as biomarkers of tumor response and overall survival 
following sunitinib treatment in gastrointestinal stromal tumors 
(GIST). Solid lines indicate relationships included in the final models, 
and dashed lines indicate relationships investigated but not included 
in the final models. 

diameters of target lesions, SLD), and overall survival from 
four clinical studies, which comprised a total of 303 patients 
with imatinib-resistant GIST, were available for PKPD analy- 
sis (Table 1).The patients had received sunitinib (25-75 mg 
orally) and/or placebo in a 4/2, 2/2, 2/1 (weeks on/weeks off), 
or continuous treatment schedule. 

Biomarker models 

The plasma concentrations of the biomarkers VEGF, 
sVEGFR-2, and sVEGFR-3 changed in a cyclic manner in 
response to therapy, returning to near-baseline levels during 
off-treatment periods, whereas the levels of sKIT continuously 
decreased over time. The biomarker time courses (BM(t)) 
were adequately described by indirect response models with 
sigmoid / max (VEGF, sVEGFR-2) or / max (sKIT, sVEGFR-3) 
drug effect relationships (Figure 2, Table 2). The maximum 
inhibitory effect (/ m ) was not significantly different from 1 
and was fixed to 1 in the subsequent analysis, assuming that 
maximum receptor inhibition was possible to achieve. Suni- 
tinib acted to decrease the zero-order production rate (KJ 
of sVEGFR-2, sVEGFR-3, and sKIT (Eq. 1 ) and to inhibit the 
degradation (K o J of VEGF (Eq. 2). The observed increases 
from baseline (biomarker at time 0, BM 0 ) in VEGF (change 
in objective function value or AOFV = 52.4) and sKIT (AOFV 
= 31 .8) concentrations over time in untreated patients were 
described by a linear disease progression model dependent 
on time (DP(t)) with a common parameter (DP s|ope ; Eq. 3). 
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The final model was simplified to include a common drug 
potency parameter (the daily sunitinib AUC resulting in half 
of the maximum drug effect, IC 50 ) for all the four biomarkers 
that were found to be highly correlated for VEGF, sVEGFR-2, 
and sVEGFR-3 (75-92%). The magnitude of interindividual 
variability in mean residence time (MRT = MK o J for VEGF, 
sVEGFR-2, and sVEGFR-3 was quantified using a common 
variability term due to the high correlation (Table 2). Data on 
sVEGFR-3 was not available for two of the studies, in total 
69 patients (Table 1). sVEGFR-2 and sVEGFR-3 were highly 
correlated for IC 50 (92%) and sVEGFR-2, sVEGFR-3, and 
VEGF were 100% correlated in terms of mean residence 
time. An individual's VEGF and sVEGF-2 levels therefore 
provided information on sVEGFR-3 responses if data were 
missing, i.e., typical values of sVEGFR-3 were not used as 
predictors for those individuals with missing data. The time 
course of sKIT was different compared with that of the other 
evaluated biomarkers with a longer mean residence time (1 01 
days vs. 3.8-23 days), and the variability in IC 50 was higher. 
An adequate predictive performance of the biomarker model 
is illustrated by a prediction-corrected visual predictive check 
(VPC) in Figure 2. A slight underprediction at later time points 
is seen for sKIT in placebo-treated patients, which may be 
due to drop out of patients with progressing disease. 

Model for tumor growth inhibition 

The data on longitudinal tumor size following placebo and 
sunitinib treatment were characterized by a tumor growth inhi- 
bition model, 12 as described by Eq. 4. The model accounts for 
the underlying tumor growth dynamics (K G ), sunitinib-expo- 
sure-driven tumor shrinkage (K ), and resistance devel- 
opment/tumor regrowth (ft(t)). In this analysis, K BM , which 
represents the tumor size reduction rate constant related to 
biomarker response, was also evaluated as predictor in the 
model. The observed tumor size at baseline (V 0 ) was incor- 
porated as a covariate acknowledging residual error (e) in the 
measurement. 19 



dy 



= K s • y(t) - (K BM ■ BM(f ) + K Druq ■ AUC) ■ R(t) ■ y(t) 



dt 

R(t) = e- 
Y(0) = Y B 



(4) 



Model-predicted sunitinib daily AUC and relative changes 
over time of sKIT and sVEGFR-3 were found to be the best 
descriptors of the change in tumor size (Table 3). In a univari- 
ate analysis, the predicted relative change in sKIT over time 
described the longitudinal tumor size data statistically sig- 
nificantly better than dose, sunitinib daily AUC, or any of the 
other biomarkers. The univariate model improved significantly 
when AUC and predicted relative change in sVEGFR-3 were 
also incorporated as multivariate predictors. The high rela- 
tive standard error for K sm (89%, Table 3) was investigated 
by log-likelihood profiling to ensure that the 95% confidence 
interval did not include zero (0.0013-0.11). 

Dropout model 

A dropout model was developed to enable prospective simu- 
lations of tumor response over time because dropout was not 
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Table 1 Data summary of the analyzed studies 


Study number 


Study 1004 


Study 1047 


Study 1045 


Study 013 


Reference 


Demetri era/. 19 


George era/. 20 


Shiaro ef a/. 21 


Maki era/. 22 


W 


202 active; initially 
47 placebo 


13 


36 


52 


Study design 


Double-blind, randomized, 
placebo-controlled, 
phase III 


Nonrandomized, evaluating 
continuous treatment regimen, 
phase II 


Nonrandomized, dose-escalating 
study in Japanese patients, 
phase l/l I 


Nonrandomized, dose-escalating 
study, phase l/l I 


Dosing schedule 


0, 50 mg q.d. 


37.5mg q.d. 


25, 50, 75 mg q.d. 


25, 50, 75 mg q.d. 


(weeks on/weeks off) 


6-week cycles (4/2) 


4-week cycles, continuous 
treatment 


6-week cycles (4/2) 


3- week cycles (2/1) 

4- week cycles (2/2) 
6-week cycles (4/2) 


Tumor assessment 
(study day) 


Cycle 1 : 0, 28 
Cycle >1 : 28 


Cycle 1 : 0, 28 
Cycles 3, 5, 7 etc: 28 


Cycle 1 : 0, 28 
Cycle >1 : 28 


2/1 : Cycle 1 : 0; Cycles 4, 8, 12, 
etc.: 14 

2/2: Cycle 1 : 0; Cycles 2, 4, 6, 
etc.: 14 

4/2: Cycle 1 : 0; Cycles 2, 4, 6, 
etc.: 28 


Biomarker assessment 
(study day) 


VEGR sVEGFR-2, 
sVEGFR-3, and sKIT 


VEGR sVEGRR-2, sVEGRR-3, 
and sKIT 


VEGR sVEGRR-2, and sKIT 


VEGR sVEGFR-2, and sKIT 




Cycle 1:0, 14,28; 
Cycle >1 : 1 , 28 


Cycle 1:0; Cycle >1: 1 


Cycle 1 : 0, 1 4, 28; Cycles 2-4: 
1, 14, 28 


2/1 : Cycle 1 : 0, 7, 14; Cycle 2: 1, 
14; Cycle >2: 1 

2/2: Cycle 1 : 0, 7, 1 4; Cycle >1 : 1 

4/2: Cycle 1: 0, 7, 14, 28; Cycle 
>1: 1 


Median tumor size at 
baseline, mm (range) 


194 (35-822) 


108 (29-191) 


166 (31-644) 


255 (55-687) 


Median survival, 
weeks (range) 


61 (4-226) 


31 (81-15) 


37 (27-48) 


39 (4-96) 



q.d., once daily; sKIT, soluble stem cell factor receptor; sVEGFR, soluble VEGF receptors; VEGR vascular endothelial growth factor. 



completely at random (those with larger tumor sizes and/or a 
poorer tumor response were more likely to drop out after the 
measurement). The final logistic regression model describ- 
ing the probability of dropping out (P) included three signifi- 
cant predictors with the estimated parameters related to time 
since start of study (0 Time ), observed SLD at dropout (0 SLD ) and 
a >20% increase in tumor size since nadir (yes/no; 9 PD ; Eq. 5, 
Table 3). P was scaled by the time interval (ATIME, typically 
6 weeks) to calculate the probability of dropout between two 
scheduled tumor assessments (P 



DR = 0, 
e 



Intercept 
DR 



Dropout' ' 

SLD+e PD PD+e Time -f 



1 + e u 



(5) 



Dropout 



1-(1-P) A 



Survival model 

All of the investigated biomarkers and tumor size measures 
were statistically significant predictors of overall survival 
when tested one by one. However, when the model-predicted 
relative change from baseline for sVEGFR-3, sVEGFR-3 REL 
(most significant predictor with AOFV of 51.6), and baseline 
tumor size were included, none of the other variables, includ- 
ing AUC, achieved statistical significance (Eq. 6). A Weibull 
model (X, hazard coefficient; a; shape factor) characterized 
the underlying baseline hazard that increased with time (a = 
1.23, Table 3, Eq. 6). Smaller decreases in sVEGFR-3 and 
larger baseline tumor size were associated with an increased 
risk of death. A clear improvement in the model fit, in terms 
of OFV, was observed when using the predicted time course 
instead of a (landmark) point estimate of change in biomarker 
levels or tumor size at 6 or 12 weeks. 



The simulation-based model evaluation comparing 
observed time course of dropout and the model predicted 
showed a reasonable description of the data (results not 
shown). A VPC of the final longitudinal tumor growth inhibi- 
tion model, taking into account dropout in the simulations, 
is presented in Figure 3 and shows good predictive per- 
formance of the tumor growth inhibition model. VPCs strati- 
fied for treatment schedule are available in Supplementary 
Figure S1 online. 



h(t) = X-a- Z"" 1 ■ gM VEGFR - 3 HEL+P 2 Timor base (g) 

Evaluation of the capability to predict survival based on 
early changes in biomarker response showed that using the 
model-predicted changes in sVEGFR-3 based on longitudi- 
nal biomarker data from only the first 6 or 12 weeks of treat- 
ment (on average, 3.1 or 4.4 measurements per individual) 
resulted in a similar fit of the survival model as when using 
the full time course of sVEGFR-3 (AOFV < 3). These results 
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Figure 2 Prediction-corrected visual predictive checks of final biomarker models based on 500 simulations (upper panels: vascular 
endothelial growth factor (VEGF), soluble VEGF receptor (sVEGFR)-2, and soluble stem cell factor receptor (sKIT); lower panels: VEGF 
placebo, sVEGFR-3, and sKIT placebo). Observed prediction-corrected biomarker data (open circles (•)), median of the observed data (solid 
line ( — )), 5th and 95th percentiles of the observed data (dashed line ( — )), and 95% confidence intervals based on corresponding percentiles 
of the simulated data (shaded areas). 



Table 2 Final biomarker model parameter estimates (relative SEs, %) 


Parameter (unit) 




VEGF 


sVEGFR-2 


sVEGFR-3 


SKIT 




Estimate 


IIV CV, % 


Estimate 


IIV CV, % 


Estimate 


IIV CV, % 


Estimate 


IIV CV, % 


BM 0 (pg/ml) 


59.8 (4.6) 


50 (8.5) 


8,660 (1.8) 


19 (6.4) 


63,900 (3.7) 


43 (12) 


39,200 (3.0) 


50 (8.1) 


MRT (days) 


3.75 (8.1) 


24 (14)" 


23.1 (4.7) 


24(14)" 


16.7 (7.2) 


24 (14)" 


101 (8.8) 


27 (25) 


Hill parameter 


3.31 (9.8) 




1.54 (10) 












IC 50 (mg-h/l) 


1 .00 (5.3) b 


50 (15) 


1.00 (5.3) b 


43 (21) 


1.00 (5.3)" 


63 (16) 


1.00 (5.3)" 


240 (14) 


DP »,„ PB </ month ) 


0.0261 (29)° 


171 (14) 










0.0261 (29) c 


172 (13) 


Residual error (%) 


44.6 (3.7) 




12.0 (8.2) 




21.9 (10) 




22.6 (8.4) 




Residual error (pg/ml) 






583 (15) 













BM 0 , biomarker level at baseline; CV, coefficient of variation; DP alo B , disease progression slope; IC 50 , daily sunitinib area under the curve resulting in half of the 
maximum drug effect; IIV, interindividual variability; MRT, mean residence time; sKIT, soluble stem cell factor receptor; sVEGFR, soluble VEGF receptor; VEGF, 
vascular endothelial growth factor. 

"Common IIV parameter. "Common IC^ parameter. "Common DP parameter. 



indicate that few sVEGFR-3 measurements in the first 12 
weeks of treatment may be sufficient to predict the probabil- 
ity of survival given the model. The VPC of the survival model 
with a Kaplan-Meier plot of the observed survival data, 
stratified by above/below median of baseline tumor size and 
sVEGFR-3 REL , overlaid with the simulated 95% confidence 
interval, shows a good predictive performance (Figure 4). 



DISCUSSION 

This PKPD analysis investigated drug exposure-biomarker- 
tumor size and survival relationships following sunitinib treat- 
ment in imatinib-resistant GIST with the aim of identifying 
relevant predictors of tumor response and overall survival. 
The soluble proteins sKIT and sVEGFR-3 were found to be 
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Table 3 Final tumor growth inhibition, dropout, and survival model parameter estimates (relative standard errors, %) 



Tumor inhibition model 






Dropout model 


Survival model 


Parameter (unit) 


Estimate 


IIV CV (%) 


Parameter (unit) 


Estimate 


Parameter (unit) 


Estimate 


K a (/week) 


0.0118 (23) 


54 (27) 


Intercept 


-3.49 (5.0) 


X (per week) 


0.00596 (49) 


Kdrus (/week/AUC) 


0.0050 (47) 


119(61) 


6 PD 


1.12 (12) 


a 


1.23 (6.9) 


K KIT (/week) 


-0.00282 (89) a 


243 (38) 


6 SLD (/mm) 


0.00105 (50) 


Pi sVEGFR-3 


3.77 (16) 


Ksvegbw < /week > 


-0.0371 (30) 






0.00707 (54) 


B OT k (/mm) 

' 2 Tumor base v ' 


0.00237 (28) 


X (/week) 


0.0217 (32) 








^cans (/ WeerS ) 


0.0017 (46) 


Residual error (%) 


12.5 (20) 








"cens 


1.27 (6.6) 



AUC, area under the concentration-time curve; CV, coefficient of variation; IIV, interindividual variability; K a , tumor growth rate constant; K mua , tumor size reduc- 
tion rate constant; K KIT , tumor size reduction rate constant; K sVEam 3 , tumor size reduction rate constant; sKIT, soluble stem cell factor receptor; sVEGFR, soluble 
VEGF receptor; VEGF, vascular endothelial growth factor; X, resistance appearance rate constant; B PD , parameter related to occurrence of disease progression; 
8 SLD , parameter related to tumor size at dropout; Q Tm , parameter related to time since start of study; X, scale factor in the Weibull probability density function; a, 
shape factor in Weibull probability density function; (5, aVEGFH . 3 , parameter relating sVEGFR-3 to the hazard; P 2Tun 



size to the hazard; X cens , scale factor in the Weibull probability density function for censoring; a ce 
ing. 

a 95% confidence interval obtained by log-likelihood profiling (0.0013-0.1 1). 



parameter relating observed baseline tumor 
shape factor in Weibull probability density function for censor- 




Figure 3 Visual predictive checks based on 500 simulations of the final longitudinal growth inhibition model for actively treated (left) and 
placebo-treated (right) patients taking into account dropout in the simulations. Observed data (open circles (•)), median of the observed data 
(solid line ( — )), 5th and 95th percentiles of the observed data (dashed line ( — )), and 95% confidence intervals based on corresponding 
percentiles of the simulated data (shaded areas). Visual predictive checks (VPCs) stratified for treatment schedule are available in 
Supplementary Figure S1 online. 



most predictive of tumor size, whereas sVEGFR-3 and base- 
line tumor size were predictive of overall survival. 

The biomarker time courses were successfully character- 
ized empirically by indirect response models where sunitinib 
inhibited the production of sVEGFR-2, sVEGFR-3, and sKIT 
and inhibited degradation of VEGF. The mechanistic rationale 
behind these biomarker modulations following sunitinib treat- 
ment is not fully understood. A proposed mechanism that is 
supported by this model is the reduction in VEGF clearance 
from blood, resulting in an increase in VEGF levels. 20 Other 
hypotheses include release of VEGF from existing pools in 
tissue or an increased production of VEGF in response to 
hypoxia. 20 GIST patients are predominantly KIT positive and 



the observed decrease in sKIT levels could be a result of 
a reduced number of viable tumor cells in treated patients. 9 
The reduced levels of circulating soluble receptor VEGFR-2 
has been demonstrated by in vitro experiments as a result 
of VEGF-mediated downregulation of VEGFR-2, 21 a putative 
mechanism also supported by this model. 

Lindauer era/. 22 described an indirect response model with 
inhibition of K m for sVEGFR-2 response in healthy volunteers 
treated with 50 mg of sunitinib for five consecutive days. In 
the current analysis, the baseline levels of sVEGFR-2 appear 
to be similar in healthy volunteers and patients, but the esti- 
mated mean residence time was four times longer, which 
may be a result of the longer treatment period. 
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Figure 4 Kaplan-Meier plot of overall survival data (blue line) and 95% prediction intervals (shaded area, 200 simulations) of the Kaplan- 
Meier plot stratified by above- and below-median baseline tumor size (upper panel) and sVEGFR-3 REL (lower panel). Median baseline tumor 
size: 195 mm, median decrease at steady state in soluble vascular endothelial growth factor receptor (sVEGFR)-3 REL : -0.32. Censoring was 
described by a Weibull model (X = 0.001 7, a = 1 .3) and applied in the simulations. 



The tumor size dynamics following sunitinib treatment were 
characterized using a previously developed tumor growth inhi- 
bition model. 12 The model has been applied to a wide range 
of clinical data. 13 ' 1416 In this example, data from placebo- 
treated patients were also described. The tumor model was 
part of a modeling framework developed to predict survival 
in phase III oncology trials based on the predicted change in 
tumor size phase II data. 12 The change in tumor size at week 
7 and the tumor size at baseline were identified as predic- 
tors of overall survival in colorectal cancer. In addition, a rela- 
tionship between survival and the model-predicted change 
in tumor size at 8 weeks and baseline tumor size has been 
established in non-small-cell lung carcinoma. 15 Our analysis 
is an extension to the previous efforts in the area. Here, the 
dynamic change in the predictors was accounted for, in con- 
trast to the models proposed by Wang 15 and Claret, 12 where 
constant value predictors were adopted. In our analysis, 
the change in tumor size was also a significant predictor of 
overall survival. However, the studied biomarkers provided a 
statistically improved description of the survival data. When 
sVEGFR-3 was related to survival, there was no additional 
improvement by incorporation of tumor size. Sunitinib primar- 
ily exhibits a cytostatic rather than a cytotoxic effect, 17 and 
change in tumor size not being an important predictor of sur- 
vival is not surprising. This may also explain why sKIT was 



the most significant descriptor of tumor size when biomarker 
response was linked to longitudinal tumor data but was not 
found to be a good predictor of overall survival. 

The finding that sVEGFR-3 response was the most 
important predictor of survival in this longitudinal model- 
based approach differs from a previous report where sKIT 
was shown to be the most important biomarker of time to 
progression and overall survival in sunitinib-treated GIST. 23 
However, the previous results were based on a traditional 
statistical analysis using discrete (landmark) time points. 
The time course of sVEGFR-2, sVEGFR-3, and VEGF are 
schedule dependent and vary widely within a treatment 
schedule. Consequently, the time point chosen to study 
the biomarker relationship will be highly influential on the 
eventual findings. In this example, we developed a modeling 
framework that allowed integration of the whole biomarker 
time course and the response. Furthermore, the use of a 
parametric time-to-event model also allowed simulations 
of the clinical outcome. The developed model enables pre- 
dictions of survival for different doses and schedules. To 
evaluate the strengths and weaknesses of these alterna- 
tive schedules, side effects should also be considered. In a 
companion article, we expand this framework to also include 
PKPD models of important side effects. 24 The established 
relationship may also enable monitoring of clinical response 
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and facilitate dose individualization because few sVEGFR-3 
measurements in the first 6 or 12 weeks of therapy appear 
to be required to predict survival given the reported model. 
Further prospective studies are, however, needed to vali- 
date the findings and to determine whether sVEGFR-3 can 
be a relevant biomarker in other sunitinib indications, such 
as renal cell carcinoma, and in treatments with other tyro- 
sine kinase inhibitors. 

Furthermore, the developed framework may be also 
applied to other molecules and indications, but the "best" 
predictors may vary and will depend on the disease and the 
mechanism of drug action. Such a framework is currently 
being constructed in renal cell carcinoma following adminis- 
tration of axitinib. 25 

In summary, an overarching modeling framework was 
developed linking exposure, biomarkers, tumor dynamics, 
and overall survival in a unified structure. sVEGFR-3 was 
found to be the most promising variable for predicting clinical 
outcome in terms of overall survival following sunitinib treat- 
ment in GIST. 

METHODS 

Patients and study design. The analysis included biomarker, 
tumor, and treatment outcome data obtained from four clini- 
cal trials in phases l-lll comprising 303 patients with ima- 
tinib-resistant malignant GIST treated with single-agent 
sunitinib 26-29 (Table 1 ). Only patients with biomarker and tumor 
data available after screening were included. Sunitinib doses 
ranged from 25 to 75 mg orally once daily and were adminis- 
tered according to one of four different treatment schedules: 
4/2, 2/2, 2/1 (weeks on/weeks off), and continuous treatment. 
The largest study (1,004) was a randomized placebo-con- 
trolled trial, where patients receiving placebo (n = 47) were 
offered sunitinib on development of disease progression as 
defined by Response Evaluation Criteria in Solid Tumors. 6 
Measurements of biomarkers and tumor sizes were collected 
for up to 1 year of treatment or until no clinical benefit. All stud- 
ies were conducted according to the Declaration of Helsinki 
and approved by the local ethics committees. Signed informed 
consent was obtained from all patients before enrollment. 

Model development. The population PKPD analysis was 
performed using the nonlinear mixed-effects modeling 
approach in NONMEM (version 7.2, gFortran version 4.5.0 
Gnu compiler collection). 30 The first-order conditional esti- 
mation method with interaction and, when appropriate, 
the Laplacian estimation method were used for parame- 
ter estimation. The R-based software Xpose (version 4) 31 
was applied for graphical diagnostics, and the PsN toolkit 
(version 3) 32 33 was used for postprocessing of results and 
execution of simulations. 

Interindividual variability was assumed to be log-normally 
distributed with a mean of zero and a variance of co 2 . Unex- 
plained residual variability was explored using additive, pro- 
portional, and combined (additive + proportional) residual 
error models or additive and combined error models on 
log-transformed scale (biomarkers). Model discrimination 
included comparison of the OFV and inspection of graphi- 
cal diagnostics. A significance level of P < 0.01 was applied 



corresponding to a decrease in OFV of at least 6.63 when an 
extra parameter was added. 

The predictive performance of the final biomarker and 
tumor model were evaluated using (prediction-corrected) 
VPCs. 34 Prediction intervals with 95% confidence intervals 
were derived from 500 simulated data sets and compared 
with the observed data. The dropout and survival model were 
evaluated by creation of Kaplan-Meier plots of observed data 
overlaid with a 95% confidence interval calculated from 200 
simulations of the studies. Imprecision (defined as standard 
errors) in model parameter estimates were obtained from a 
nonparametric bootstrap with 200 samples. 

The PKPD model-building process resulting in the final 
model linking the predictors to overall survival was developed 
in a sequential manner. The different steps are described 
below. 

Pharmacokinetics. Individual parameter estimates or, when 
no PK information was available (57 patients), typical pop- 
ulation parameters obtained from a previously developed 
PK model 35 were used to describe the concentration-time 
profiles of sunitinib (no data for the equipotent metabolite 
SU1266 was available) Information on protein binding was 
also not available. Dose, daily AUC with AUC = 0 during off- 
treatment periods, or predicted concentration-time profiles 
were evaluated as exposure measures. AUC was calculated 
as Dose/(CL/F), where CL/F is the apparent oral plasma 
clearance. 

Biomarker models. Biomarker modulations in placebo- 
treated patients due to the natural history of the disease were 
explored and characterized using linear or nonlinear disease 
progression models. Changes in biomarker concentrations 
over time were related to sunitinib exposure by indirect 
response models with linear, inhibitory E max (/ ), and sig- 
moid / max drug effect relationships. 36 Modeling was performed 
on log-transformed biomarker data. 

Models for each of the biomarkers were initially devel- 
oped separately and finally combined into a joint model to 
explore correlations between the biomarkers and to simplify 
the model structure (where appropriate). For the joint model, 
AUC was used as a measure of drug exposure because the 
complexity of the model resulted in long computational run 
times. Because the independent parameter estimates for 
disease progression (DP s|ope ) and the drug exposure result- 
ing in half of / max (IC 50 ) were similar among the biomarkers 
when estimated separately, they were evaluated as shared 
(common) parameters (see Supplementary Data). 

Tumor growth inhibition model. The relationship between the 
biomarkers and tumor size (SLD) was characterized by a 
tumor growth inhibition model. 12 Other tumor growth models 
were evaluated, including models with Gompertz equation 
and zero-order tumor growth, but they did not provide a bet- 
ter description of the data. Sunitinib dose, sunitinib daily AUC 
the individual model-predicted relative or absolute change 
from baseline over time, and the absolute value for the differ- 
ent biomarkers (BM(t)) were tested separately or in combina- 
tion to describe the response. 

Dropout model. Prospective tumor model simulations required 
the frequency and time course of patient dropout from tumor 
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measurements to be taken into account. The exact date of 
dropout was not available (unless because of death) and 
was therefore assumed to occur at the time of last available 
tumor assessment. Observed and predicted baseline SLD, 
SLD at time of dropout, difference in SLD to baseline, and 
the change in SLD at dropout since the last assessment 
were evaluated as factors affecting the probability to drop out 
over time in a logistic regression model. In addition, a >20% 
increase in SLD since nadir, time since start of study, dose, 
AUC, and study were tested for significance. In the simula- 
tion-based model evaluation, dosing records were imputed 
according to protocol until the time of the last observed tumor 
assessments in each substudy (see Supplementary Data). 

Overall survival model. A parametric time-to-event model 
was developed to explore the relationships between the 
four biomarkers, tumor size, and overall survival. The 
exponential, Weibull, log-logistic, extreme value, and Gom- 
pertz probability density functions were compared to best 
describe the observed survival data. Dose, AUC, model- 
predicted baseline levels of the biomarkers and tumor size, 
and the relative change over time in tumor size and bio- 
markers were evaluated as predictors of overall survival. 
Although there would have been potential value in using 
the full PK profile as driver for efficacy, the complexity of 
the framework and the resultant run times severely lim- 
ited the practicalities of such an approach. The predictors 
were extrapolated based on developed models until time 
of death/censoring. Furthermore, a point estimate of the 
relative change in biomarker levels, in addition to the tumor 
size at treatment week 6 or 12 (i.e., after completion of one 
or two treatment cycles), was evaluated to be related to the 
survival data. In the simulations, censoring, e.g., because 
of a short follow-up period, was described by a Weibull 
model (see Supplementary Data). 

Acknowledgments. E.K.H. was supported by the Swedish 
Cancer Society, Sweden. Parts of the computations were 
performed on resources provided by the Swedish National 
Infrastructure for Computing (SNIC) at UPPMAX. Martin Ag- 
back at UPPMAX is acknowledged for assistance concerning 
technical aspects in making NONMEM run on the UPPMAX 
resources. 

Author contributions. E.K.H. , M.A.A., P.A.M., J.F., M.O.K., 
and L.E.F. wrote the manuscript. E.K.H., M.A.A., P.A.M., 
B.E.H., J.F., M.O.K., and L.E.F. designed the research. 
E.K.H., M.A.A., P.W., J.F., M.O.K., and L.E.F. performed the 
research. E.K.H. and P.W. analyzed the data. 

Conflict of Interest. M.A.A., B.E.H., and P.A.M. are employ- 
ees of Pfizer Ltd. At the time the work was carried out, J.F. 
was an employee at Pfizer Ltd. L.E.F. and M.O.K. have acted 
as paid consultants to Pfizer Ltd. P.W. was supported by a 
research grant from Pfizer Ltd. As Deputy Editor-in-Chief for 
CPT: Pharmacometrics & Systems Pharmacology, L.E.F. 
was not involved in the review or decision process for this 
article. The other authors declared no conflict of interest. 



Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y Models describing longitudinal tumor size data 
are being increasingly applied and the tumor 
size at a specific time point is linked to survival. 
For antiangiogenic drugs, such as sunitinib, the 
time course of biomarkers may offer greater 
predictive value for clinical response. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

S This article describes the development of an 
overarching modeling framework linking drug 
exposure, biomarkers, tumor dynamics, and 
overall survival in a unified structure. The time 
course of variables was tested for their predic- 
tive capacities to identify predictors of clinical 
response. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

S The relative change in sVEGFR-3 over time and 
tumor size at start of treatment were found to be 
the most promising variables for predicting clini- 
cal outcome, with decreased sVEGFR-3 levels 
and smaller baseline tumor size being predic- 
tive of longer overall survival. 



HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S The developed framework may be used to per- 
sonalize anticancer treatment by facilitating 
dose individualization, monitoring of clinical 
response, and the identification of patients most 
likely to benefit from treatment. 
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